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role in the determination of the j5 s 
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Abstract. This talk consists of two parts. We first present a light-cone QCD sum rule computation 
of the B s — > /o(980) form factors which are necessary inputs in semileptonic and nonleptonic 
B s decays into /o(980). Then we analyze nonleptonic B s decays into a charmonium state and 
a light meson, which are potentially useful to access the B s -B s mixing phase j3 s . We explore 
the experimental feasibility of measuring these various channels, paying attention to different 
determinations of j8 s in view of the hints of new physics recently emerged in the B s sector. 
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INTRODUCTION 

The detailed study of CP violation is a powerful and rigorous tool in the discrimination 
between the Standard Model (SM) and alternative scenarios. For instance the analysis of 
the B s unitarity triangle of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements: 
VusV^ + V cs V* b + VtsV* b = provides an important test of the SM description of CP 

violation. One of its angles, defined as fi s = Arg [ — y^fi ] , is half of the phase in 

the B s -B s mixing, and is expected to be tiny in the SM: f5 s ~ 0.017 rad. The current 
measurements, by the CDF and D0 collaborations at Tevatron based on the angular 
analysis of the time-dependent differential decay width in the process B s — > J /yty [1], 
indicate larger values and the averaged results are consistent with SM only at 2.2 o level: 
^v/v* = _ 2 p s = _ .77±0.29 or qJ/y* = _ 2 p s = _2.36±Jg [2]. Although the recent 
result by the CDF: /3, G [0.0,0.5] U [1.1, 1.5] (at 68% confidence level) [3] has a smaller 
deviation from the SM, the uncertainties are still large and the precise measurement of 
/3 S is a priority for the forthcoming experiments. Towards this direction the nonleptonic 
B s decays are certainly of prime importance. 

In this work we first compute the B s — > /o(980) form factors using the light-cone QCD 
sum rule (LCSR)[4]. These results will be useful in the analysis of semileptonic and 
nonleptonic B s — > fo decays. Subsequently we investigate the B s decay modes induced 
by the transition b — >■ ccs, namely B s — >■ M c g + L, where M cc - is an s-wave or p-wave 
charmonium state and L is a light scalar, pseudoscalar or vector meson, /o(980), r}, r\' ', 
[5]. In particular, we exploit the generalized factorization approach to calculate their 
branching fractions in the SM in order to understand which of these modes are better 
suitable to determine /3 S . 



B s ->> f FORM FACTORS IN LCSR 

Hereafter we will use /o to denote /o(980) meson for simplicity. The parametrization of 
matrix elements involved in B s — > /o transitions is expressed in terms of the form factors 
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(MPfM\Bs(PB s )) = -i{Fi(q 2 )[P,- *' fQ q,} + F (q 2 ) * fo q,}, 
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where P = p Bs + p fo , q = p Bs - p fo , and J 5 ^ = sy^y 5 b, J 5 * = so^ v y 5 q v b . To compute 
such form factors in the LCSR [6] we consider the correlation function: 

n(p fQ ,q)=ij d 4 xe^(M Pfo )\T{j ri (x)Jr 2 (0)}\0} (2) 

with jr x being one of the currents in the definition of the B s — >■ /o form factors: jr l = 
for Fi and Fq, and jr x = Jf? for Ft. The matrix element of jr 2 = biy^,s between the 

2 

vacuum and B s defines the B s decay constant f Bs : (B s (p Bs )\biy 5 s\0) = j^^/b,- 

The LCSR method consists in evaluating the correlation function in Eq. (2) both at 

the hadron level and at the quark level. Equating the two representations allows us to 

obtain a set of sum rules suitable to derive the form factors. 

The hadronic representation of the correlation function in Eq. (2) can be written as 

the sum of the contribution of the B s state and that of the higher resonances and the 

continuum of states h: 

11 {Pfo,q) = 2 — 7 — ; — ^2 / ds — ? — ; — \2> ( 3 ) 

™ Bs -(pfo+q) A) s-(pfo+q) 

where higher resonances and the continuum of states are described in terms of the 
spectral function p h (s,q 2 ), contributing above a threshold sq. 

The correlation function can be evaluated in QCD with the expression 



Im nQ CD ( 5 ,4 2 ) 

K+m,) 2 " s-(p fo +q) 2 



(4) 



Expanding the T-product in Eq. (2) on the light-cone, we obtain a series of oper- 
ators, ordered by increasing twist, the matrix elements of which between the vac- 
uum and the fo are written in terms of /o light-cone distribution amplitudes (LCDA). 
Since the hadronic spectral function p h in (3) is unknown, we use the global quark- 
hadron duality to identify p h with pQ CD = ilmnQ CD when integrated above sq so that 

f£ ds s ffr*+ q )2 = hG ds ^0jfT7jp- - Usin § the quark-hadron duality, together with 

the equality Tl n (pf Q ,q) = Yl^ CD (pf Q1 q) and performing a Borel transformation of the 
two representations, we obtain a generic sum rule for the form factors 

_ _ _ m B s 1 fSO _ s 

(M Pfo )\j ri \B s (p Bs ))(B s (p Bs )\jr 2 \0)e ^ = - / ds e IP lmTl^ CD ( s ,q 2 ), (5) 
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TABLE 1. B s -> f form factors in the LCSR. 
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where Q 2 = —q 2 , pb s = Pf +q an d M 2 is the Borel parameter. The Borel transformation 
will improve the convergence of the series in n^ CD and for suitable values of M 2 
enhances the contribution of the low lying states to I1 H . Eq. (5) allows us to derive 
the sum rules for F\, Fq and Fj, choosing jr l = or jy x = J^ T ■ 

We refer to Ref. [4] for numerical values of the input parameters as well as for the 
final expressions of the form factors obtained from (5). The sq is supposed to be around 
the mass squared of the first radial excitation of B s and is fixed as so = (34 ± 2) GeV 2 . 
As for the Borel parameter, the result is obtained requiring stability against variations of 
M 2 . In Fig. 1 we show the dependence of F\ (q 2 = 0) and Fr(q 2 = 0) on M 2 and we find 
the stabilities when M 2 > 6 GeV 2 , and thus we choose M 2 = (8 ±2) GeV 2 . 

To describe the form factors in the whole kinematically accessible q 2 region, we use 

the parameterization F^q 2 ) = 2 2 , m i vz > z G f 1 ' ' T ^ the P aramQtQTS 

a\ and bi are obtained through fitting the form factors computed numerically in the large 
recoil region. Our results are collected in Table 1, where uncertainties in the results are 
due to the input parameters, sq and M 2 . 



B s M C cL DECAYS 



The effective hamiltonian responsible for decays induced by the b — >■ ccs transition is: 

^ eff = -^{v cb v c * s [c 1 (ii)o l +c 2 (ii)02\-v tb v; s [ £ Ci(ji)Oi(ji)]}, (6) 

where Gf is the Fermi constant, Oi are the four-quark or magnetic-moment operators 
and C, are Wilson coefficients. With the assumption of the CKM unitarity and the neglect 




of the tiny V ub V u * s , we have V tb V* = -V cb V* s . 

The simplest approach to compute the matrix element of the effective four-quark 
or magnetic-moment operators is the naive factorization approach. In this approach 
neglecting the magnetic moment operators whose contributions are suppressed by a s , 
the B a — y M C cL amplitude reads ((a being a light flavour index)) 

srf{B a M C ,L) = ^L Vcb V*af(fi) {M cd \cf{\ - y 5 )c \0)(L\sy^l - y 5 )b\B a ) , (7) 

where a| ff (ju) is a combination of the Wilson coefficients: «| ff (/i) = «2(m) +a^(jl) + 
as(n) and and a.2 = C2 + C\/N c , a 3 = C^-\-CA/N c + \e c (C 9 + Ci /N c ) and a 5 = C 5 + 
C(,/N c + \e c (C-j + C^/Nc) with e c = 2/3 and N c = 3. This factorization approach allows 
us to express the decay amplitudes in terms of the heavy-to-light form factors and the 
decay constant of the emitted meson. Unfortunately, one severe drawback is that naive 
factorization badly reproduces several branching ratios for which experimental data are 
available. In particular, the b — y ccs induced modes under scrutiny are color suppressed, 
and the predictions of naive factorization will typically undershoot the data. The most 
striking discrepancy is for the Bj decay modes with Xc-o in the final state, which have a 
sizeable rate but their decay amplitude vanishes in the factorization approach [7]. 

Several modifications to the naive factorization ansatz have been proposed and in 
particular in this work we will explore one possibility by treating the Wilson coefficients 
as effective parameters to be determined from experiment. In principle, it implies that 
such coefficients are channel dependent. However for channels related by invoking 
flavour symmetries, universal values for the coefficients can be assumed. In our case, 
this generalized factorization approach consists in considering the quantity af f in (7) as 
a process dependent parameter to be fixed from experiment. In particular, if one assumes 
the flavor St/ (3) symmetry, B (B u or Bj) decays can be related to analogous B s decays, 
so that experimental data on B decays provide a prediction for B s related ones. 

Our strategy is to exploit the existing experimental data for B decay modes to deter- 
mine an effective parameter af 1 and, assuming SU (3)p symmetry, to use these values to 
predict the flavour related B s decays. In the case of modes with % c \ in the final state, we 
will determine the combination f Xcl a| ff since sizable uncertainties may be introduced to 
the Wilson coefficient but will cancel in the predictions of branching ratios. In this pro- 
cedure, we use two sets of form factors: the one obtained using sum rules based on the 
short-distance expansion [9], and the set in [10] based on the light-cone expansion. In 
the case of B s — y <j> and B s — y /o(980) we use form factors determined by LCSR [10, 4]. 
B s —y 7]W form factors are related to the analogous B —y K form factors and the mixing 
angle between r\ and r\' in the flavor basis [11] can be fixed to the value measured by 
the KLOE Collaboration: 6 = (41.5 ± 03 stat ± 0.7. v> , st ± 0.6 th ) ° [12], which is also sup- 
ported by a QCD sum rule analysis of the radiative — y TjWy modes [13]. The Wilson 
coefficients for B s — y M c -c/o(980) are obtained using the effective value determined from 
B M C cK. Most of other numerical inputs are taken from the particle data group and 
we refer to Ref. [5] for more details. 

The predictions for branching ratios of B s decays are given in Tables 2 and 3. In Table 
2 the available experimental data [8, 14, 15] are also reported, with a satisfactory agree- 
ment with the predictions. In theoretical predictions we have included the uncertainty 



TABLE 2. Branching ratios (in units of 10 4 ) of B s — > M c gL using the form factors in [9] 
(CDSS) and in [10] (BZ). Experimental results are taken from PDG [8], except for B s -> 
7/i/A77 (77') measured by Belle Collaboration [14] and the bound for B, -^J/wfo is fromRef. [15]. 
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TABLE 3. Branching ratios of B s decays into p-wave charmonia (unit: 


10" 4 ). 
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on the form factors at q 2 = and on the experimental branching ratios, but in the case 
of the modes involving r\ or r\ ' the uncertainty on the form factors is not included since 
the dependence on the form factors will cancel when the branching ratios of B — > J / yK 
decays are related to the corresponding B s decays. 

As appears from Tables 2 and 3 all the considered modes have sizable branching frac- 
tions which are large enough to make them promising candidates for the measurement of 
fi s . The modes involving r\ , 7]', fo present, with respect to the golden mode B s — >■ 
the advantage that the final state is a CP eigenstate, not requiring angular analysis. How- 
ever, channels with r\ and r\ ' can be useful only after a number of events have been 
accumulated, since at least two photons are required for the reconstruction. 

As discussed in [16, 4, 17], B s — > J/tyfo has appealing features since, compared 
with the tj( ), the fo can be easily identified in the k + tz~ final state with a large BR: 
3§{ fo — > Tl + Tl~) = (50+g)% [18], so that this channel can likely be accessed. At present, 
the Belle Collaboration has recently provided the following upper limit [15]: 

@{B S -> J /Wo) x ^(/o -> < 1-63 x 1(T 4 (8) 

marginally in accordance with our prediction. 

Let us come to B s decays to p-wave charmonia. Among these decays, the only one 
with non vanishing amplitude in the factorization assumption is that with % c \ in the final 



state. In the other cases, i.e. modes involving Xc-o.2 or h c , which we show in Table 3, 
results are obtained determining the decay amplitudes from the B decay data by making 
use of the SU(3) symmetry. In this case, the differences between the B and B s decays 
arise from the phase space and lifetimes of the heavy mesons. As for the mechanism 
inducing such processes, one possibility is that rescattering may be responsible of their 
observed branching fractions, as proposed in Ref.[7]. Among these channels, B s — > XcO<l> 
is of prime interest and promising for both hadron colliders and B factories. 

CONCLUSION 

Recent results in the B s sector strongly require theoretical efforts to shed light on which 
are the most promising decay modes to unreveal new physics. In this work we have 
analyzed channels induced by the b — > ccs transition. Modes with a charmonium state 
plus r\, r\' , /o(980) are the most promising, being CP eigenstates not requiring an 
angular analysis. In particular, the case of /o is particularly suitable in view of its easier 
reconstruction in the subsequent decay to k + k~ . As a preliminary step we have used 
the light-cone sum rules to compute the B s — > /o(980) form factors which are necessary 
inputs in the analysis of B s decays. 
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